A densest compact planar packing 
with two sizes of discs 

Tom Kennedy 
Department of Mathematics 
University of Arizona 
Tucson, AZ 85721 
email: tgk@math.arizona.edu 

February 1, 2008 

Abstract 

We consider packings of the plane using discs of radius 1 and r = 0.545151 • • •. 
The value of r admits compact packings in which each hole in the packing is formed 
by three discs which are tangent to each other. We prove that the largest density 
possible is that of the compact packing shown in figure ^ 
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1 Introduction 



We consider the following packing question in two dimensions. Fix a number r < 1. Using 
discs of radius 1 and r, what is the densest packing of the plane? We do not impose any 
constraint on the relative number of discs of the two sizes. 

It was proved long ago that the densest packing of the plane using discs of equal radii 
is to put the centers of the discs on a triangular lattice. [201 1^- (The Voronoi cells of this 
packing are hexagons, and this packing is often referred to as hexagonal or honeycomb.) 
The density of this triangular packing is tx j^Vl. L. Fejes Toth observed that if r is slightly 
less than 1, then one cannot do any better than this packing density [S]. The interval 
in which it has been proved that the highest packing density is 7r/v^T2 was increased to 
[0.906 ■ ■ ■ , 1] by Florian [TT]. Using an an idea of Boroczky [Bo], Blind [21 E] and G. Fejes 
Toth 5j independently extended it to [0.742 ■ ■ ■ , 1]. 

For smaller values of r there is a rich variety of packings with densities greater than 
7r/vT2. (A survey of the best known packings as a function of r may be found at |17;.) 
However, the densest packing has been rigorously established only for six particular values 
of r. All six of these values of r allow compact packings. A packing is said to be compact 
if each disc is surrounded by a ring of discs, all of which are tangent to the disc at the 
center. Furthermore, each disc in the ring is tangent to the two discs in the ring which 
are adjacent to it in the cyclic order. Heppes has proved that for six values of r which 
allow compact packings, the largest density is attained by a particular compact packing. 
It has been shown that there are only nine values of r that admit compact packings 

In this paper we consider one of the values of r which admits compact packings but 
for which it has not been shown that a compact packing attains the largest density. The 
value is r = 0.545151042 ■ ■ ■. The exact r is a root of 

- 8r^ - 44r^ - 232r^ - 482r^ - 24r^ + SSSr^ - 120r + 9 = (1) 

(This equation is derived in appendix O) In the remainder of this paper r will denote 
this particular radius. In this paper we prove that the largest packing density possible 
using discs of radius 1 and r is that attained by the compact packing in figure It has 
a packing density of 5 = 0.911627478 ■ ■ ■. 

The strategy of our proof comes from a technique in classical statistical mechanics 
known as "m-potentials" ^Hl- Our strategy is similar to Heppes "cell balancing" JH]- In 
statistical mechanics m-potentials were introduced to deal with frustrated spin systems. 
In spin systems the Hamiltonian or energy function is typically a sum over translates of 
a local energy function. A ground state of such a system is a configuration of spins which 
minimizes the total energy function. One can also ask what is the minimum of a single 
local energy function. This may be less than the value of the local energy function in the 
ground state. When this happens the system is said to be frustrated. What one would like 
to do locally to minimize the energy cannot be done globally to simultaneously minimize 
all the local energy functions. 
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Figure 1: A densest packing with r = 0.545 ■ ■ ■. The quadrilateral shows a unit cell for 
the packing. 

In classical spin systems an m-potential is a local function on the spin configurations 
with the property that when it is summed over all translates the result is just the zero 
function. Thus one can add this local function to the local energy function and obtain a 
total energy function with the same ground states. In some problems by carefully chosing 
the m-potential one can obtain a system which is not frustrated. (We should note that 
an important class of frustrated spin systems comes from disordered systems such as spin 
glasses, but m-potentials have not been useful in this context) 

The disc packing problem is also frustrated. Given just three discs the densest packing 
is to have them all touch one another. For most values of r it is not possible to find a 
packing in which the packing locally always consists of three discs touching each other. 
Even when it is possible (as it is for our value of r), different choices for the three discs 
that are tangent will give different densities for the triangle formed, and it is not possible 
to find a packing that only uses the densest triangle. By analogy with the m-potentials, 
we will introduce a function on packings that is a sum over all triangles of a function of 
three discs. We will refer to it as a "localizing potential" since its purpose is to reduce 
the global problem of finding the best packing to a local problem involving only three 
discs. Note that for the spin problem, if we are only interested in the ground states then 
it suffices that the sum of the local m-potentials be non-negative. Likewise, in the packing 
problem it suffices that the sum of the localizing potentials be non-negative. 
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In this paper we only allow discs of two sizes. An interesting and presumably more 
complicated question is what is the largest packing density if we allow discs with any 
radius in [r, 1]. Another interesting question is what is the densest packing if we add the 
constraint that the ratio of the number of one type of disc to the other must converge 
to a given value in the limit of packing the entire plane. This question was studied 
non- rigorously in ^3]. 

In the next section we explain the method of localizing potentials in detail. Our 
localizing potential is the sum of two parts, a "vertex localizing potential" and an "edge 
localizing potential." In section IHl we define the vertex localizing potentials used in the 
proof of our result and show that the sum of our vertex localizing potentials is non- 
negative. We define the edge localizing potentials and show their sum is non-negative in 
section m Sections El and IHl are devoted to studying the local problem involving just three 
discs that comes from adding our localizing potential to the original packing problem. 



2 Localizing potentials 

Consider the centers of the discs in a large packing. The Delaunay decomposition gives a 
triangulation in which the vertices of the triangles are the centers of the discs. We denote 
the triangles by Tj. Let A{T) be the area of triangle T. Given a triangle T, let 0i) 02 
be the angles in T and rj the radii of the discs at the corresponding vertices. We define 



DiT) = f2l<j^^r^ (2) 



2 

i=0 

As long as the triangle is not too "flat", D{T) is the area of the intersection of the triangle 
with the three discs. Even when the triangulation contains triangles for which this is not 
true, the sum ^j-D(7i) will be (up to boundary effects) the total area of the discs in 
the packing since the sum of the angles around a vertex is always 2tt. Up to boundary 
effects the packing density is 



So we want to prove this ratio is no greater than 5. This is equivalent to 



(3) 



Y,i^A{T;) - D{T,)) > (4) 

i 

We define 

E{T) = 5A{T) - D{T) (5) 

Heppes defines the "surplus area" to be A{T) — D(T)/6. We refer to the quantity E(T) 
as the "excess" of the triangle. If it were nonnegative for every triangle we would be done. 
This happens to be true when one considers packings with discs of a single radius, but 
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it is not true in our problem. In an optimal packing, triangles have both positive and 
negative excess, but the sum over all triangles of the excess is zero. We must prove that 
for any packing the sum over the triangles of the excess is non-negative. 

We want to define a function F{T) on triangles with the following two properties. 
First, for any Delaunay decomposition we require 

^F(T,)>0 (6) 

i 

Second, for any triangle that can occur in a Delaunay decomposition, we require 

E{T) - F{T) > (7) 
If we can do this, then we are done : 

j2m)>J2(^m)-Fm>o (s) 

i i 

We refer to F as a localizing potential since it reduces proving the global inequality (jHl) to 
proving the local one (j?!). We will prove the following theorem. 

Theorem 1 For the particular value of r we are considering (r 0.545151j, there is a 
localizing potential F{T) which satisfies inequalities ^ and ^ with 6 equal to the density 
of the packing shown in figureU\ Thus the density of a packing consisting of discs of radius 
1 and r is at most 6. 

We parameterize triangles by their edge lengths and the radii of the discs at these 
vertices. We label the vertices i = 0, 1,2 and we denote the length of the edge opposite 
vertex i by Xi. The radius of the disc at vertex i is rj. The excess is then written as 
E{xo, xi, X2, To, Ti, r2), and the localizing potential as F{xo, xi,X2, tq, ri, r2). F will be the 
sum over the three vertices of the triangle of a vertex localizing potential plus the sum 
over the three edges of an edge localizing potential. The vertex potential is based on the 
constraint that the sum of the angles around a vertex is 27r. The edge potential is based 
on a constraint involving the signed distance from the edge to the center of the circle 
which circumscribes the triangle. 

Consider the vertex potential v for vertex 0. It depends on xq, xi,X2 only through the 
angle of the triangle at vertex 0, which we denote by (pQ. So we write the vertex potential 
as f (00, To, ri, r2). The vertex potentials for the other two vertices are w ri, tq, 
for vertex 1 and v{(f)2,r2,ro,ri) for vertex 2. We will always take v{4>o,ro,ri,r2) to be 
symmetric under the interchange of its last two arguments. 

Now consider the edge localizing potential for the edge opposite vertex 0. It will be 
a function of the edge lengths Xq, Xi, X2 and the radii of the discs at the endpoints of the 
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edge, i.e., ri, r2. We write it as e(xo, xi, X2, ri, r2). It is symmetric under the simultaneous 
interchange of ri with r2 and Xi with X2. The total localizing potential for our triangle is 

F{xo, Xi,X2, ro, ri, r2) = f (0o, ^0, n, ^^2) + n, ?^o, ^2) + v{4>2, r2, tq, ri) 

+e(xo, xi, X2, ri, ra) + e(xi, Xq, X2, ro, r2) + e(x2, Xq, Xi, tq, n) (9) 

Note that each edge in the packing has two edge locahzing potentials associated with 
it, while each vertex has n vertex localizing potentials associated with it where n is the 
number of triangles in the packing that contain the vertex. Thus the sum of F{T) over 
all triangles can be written as the sum of the following two sums. The first is the sum 
over edges of the sum of the two edge localizing potentials associated with the edge. The 
second is the sum over vertices of the sum of the n vertex localizing potentials associated 
with the vertex, {n is vertex dependent.) Thus to prove (jH}, it suffices to prove the 
following two conditions. First, for every disc center in the packing we have 

n 

J2^i^i^r,ri,ri+i)>0 (10) 

i=l 

where n > 3 is the number of triangles that have the disc center as a vertex, r is the 
radius of the disc, and ^i, r2, ■ ■ ■ , r„ are the radii of the discs that surround it. These are 
ordered in the natural way, so that one triangle has discs of radii r, ri,r2, the next has 
radii of r, r2, and so on. r„+i is defined to be ri. 6'i, 6*2, ■ ■ ■ , 6'„ are the angles of these 
triangles at the vertex. Second, for every edge in the packing we have 

e(xo, xi, X2, ri, r2) + e(xo, x'^, 4, n, ^^2) > (11) 

Here the length of the edge is Xq, and ri, r2 are the radii of the two discs at its endpoints. 
xi,X2 are the lengths of the other two edges in one triangle and the lengths 

of the other two edges in the other triangle. The edge potential e(xo, xi, X2, ri, r2) will 
depend on the signed distance from the "center" of the triangle to the edge opposite 
vertex 0. The center of the triangle is the point equidistant from the three vertices. The 
signed distance is positive when the center lies on the same side of the edge as the vertex 
opposite the edge. For two triangles in a Delaunay decomposition that share an edge, 
there is a constraint on the two signed distances to this common edge that will be the 
basis for proving (fTT|) . See section HJ 

We end this section with a discussion of why we use the Delaunay triangulation rather 
than the FM triangulation 9^ that was used by Heppes in his proof of the optimality of 
six other compact packings. For most triangles the excess is positive. As we will see in the 
next section, it can be negative for triangles that are close to some of the triangles that 
appear in the compact packing of figure Q The only other triangles with negative excess 
are relatively fiat triangles, i.e., triangles with a large obtuse angle. The FM triangulation 
has the nice property that each disc is covered by the triangles that have a vertex at the 
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center of the disc, and so triangles cannot be too flat. However, an FM triangulation can 
still contain triangles which are flat enough that their excess is slightly negative. The 
Delaunay triangulation can contain triangles that are even flatter and so have an even 
more negative excess. But in the Delaunay triangulation it is possible to define an edge 
localizing potential that takes care of these flat triangles. This is explained in section El 



3 The vertex localizing potentials 

We begin the proof of the theorem by considering what the values of the localizing poten- 
tials must be for the triangles that appear in our densest packing. There are four triangles 
that appear in this packing. We will refer to the triangle that has two large discs and one 
small disc as the alpha triangle. It has one side of length 2 and two of length 1 + r. We 
denote the angle opposite the side of length 2 by a and the other two angles by a! . For 
this triangle the excess is negative, 

^ -0.0022743457 (12) 

We will refer to the triangle that has one large disc and two small discs as the beta 
triangle. It has one side of length 2r and two of length 1 + r. The angle opposite the side 
of length 2r will be called (3, and the other two are (3'. This triangle also has negative 
excess. 

Ef3 ^ -0.0017217279 (13) 

The triangle with three small discs will be called the small equilateral triangle. Its 
sides are of length 2r, and its excess is positive. We denote it by 

Es ^ 0.0024336170 (14) 

The fourth triangle has three large discs and will be called the large equilateral triangle. 
Its sides are of length 2, and its excess is also positive. We denote it by 

El ^ 0.0081887688 (15) 

The localizing potential must be defined so that it equals the excess for each of the 
four triangles. The edge localizing potential is zero for all four of the triangles. So we 
have the following four conditions. 

v{a,r,l,l) + 2v{a',l,r,l) = E^ (16) 

v{(3,l,r,r) + 2v{(3',r,r,l) = Ef, (17) 
3v{7T/3,r,r,r) = Es (18) 
3v{7t/3,1,1,1) = El (19) 
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The localizing potential condition ()10|) for a small disc requires 

2v{a, r, 1, 1) + 2v{(3', r, 1, r) + f (tt/S, r, r, r) = (20) 

and for a large disc it requires 

v{(3, 1, r, r) + Av{a', 1, r, 1) + 2v{tt/3, 1, 1, 1) = (21) 

We have found six conditions above, but these six conditions are not linearly indepen- 
dent. In the packing shown in figure ^ the unit cell has 6 alpha triangles, 3 beta triangles, 
1 small equilateral triangle and 2 large equilateral triangles. The sum of the excesses of 
the triangles in a unit cell must be zero, so 

6Ea + 3E(3 + Es + 2El = (22) 

This implies that the above six conditions are equivalent to five linearly independent 
conditions. 

We solve these equations by first introducing two parameters x and y. We set 



f (a, r, 1, 1) = xEa 

!;(«', l,r,l) =i(l-a;)E„ 
v{f3,l,r,r) =yEp 

T;(/3',r,r,l) =\{l-y)Ep 
v{Tc/3,r,r,r) = 

t;(7r/3, 1,1,1) =IEl (23) 

The conditions ()20|) and (PT|) are equivalent conditions on x and y. They give 

x2E^ - yE(s = 2E^ + '^El = -E^ - ^Es (24) 

We take x = 0, and then y is determined by the above equation. We will explain the 
motivation for this choice at the end of this section. 

We let 00(^^0, ''^1, ''^2) denote the angle 0o in the triangle with discs of radius ro,ri,r2. 
So it equals a,a',f3,(3' or 7r/3. Then the vertex localizing potential is 

v{(f),ro,ri,r2) = t^(0o('^o, n, ^^2), ^o, ri, r2) +m|0 - 00(^^0, n, ^^2)! (25) 

where f (00(^^0, ''^1, ''^2), ''"0, ''^1, ''"2) is given by the above equations. So equations (j2Sl) are 
minima of the localizing potential, and the potential increases with slope ±m as we move 
away from these minima. Note that we use the same slope, m, for all the localizing 
potentials. We will take m = 0.12. 
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Next we prove that the vertex locahzing potential we have defined satisfies ()10j). We 
consider a disc D and must show 



^t;(^i,p,ri,ri+i) > (26) 

2 = 1 

Here n > 3 is the number of triangles with a vertex at the center oi D. p is the radius of 
D, and ri, r2, ■ ■ ■ , r„ are the radii of the discs that surround D. The angles 6'i, 6*2, ■ ■ ■ , 6*^ 
are the angles of these triangles at Of course, the sum of the angles around a vertex 
is 2tt. So 

n 

J]^2 = 27r (27) 
1=1 

For p = r, let rirrr be the number of i with rj = rj+i = r. Let Urii be the number of i 
with = Tj+i = 1. And let n^ir be the number of i with one of r^, rj+i equal to r and one 
equal to 1. (Keep in mind the symmetry f (0, Tq, ri, r2) = f (0, Tq, r2, ri).) Then to prove 
inequality for p = r it suffices to prove 



Eg 1 — y TT 

rirrr 1" '^rrl Ef} + TlrllXEa + m\2Tl — Urrr UrrlP' — nrUC(\ > (28) 

3 2 3 

And to prove the inequality for p = 1 it suffices to show 

riirryEp + nirJ—^Ea + '^ni"^ + "^|27r - ni„(3 - niria - '^iii^l > (29) 

where ni„i ni^-i and rtm are defined in the obvious way. Each choice of n^rr, n^n gives 
a lower bound on m, as does each choice of nirr, 'n.iri, "^.m. Note that for n„j. = = 
2, ririi = 2, the left side of is zero, and for riirr = 1, "n-iri = 4, nm = 2, the left side of 
fl29|l is zero. These are the cases which occur in the optimal packing. The largest lower 
bound on m that we find from the other cases is when nirr = 4, niri = 4, nm = 0. This 
case implies 

m> 0.1185912 (30) 

We take 

m = 0.12 (31) 

The above computation is the motivation for the choice of a; = 0. Other choices force m 
to be larger. 

Finally, we make a modification to definition We define f (0, tq, ri, r2) to be 

given by the above definition provided the value is less than 0.1. Otherwise we define it 
to be 0.1. Proving (jlOp for this modified function is easy. The most negative value of 
w(0,ro,ri,r2) is E^ ^ -0.0022743457. So if one or more of the v in dTUD is 0.1, then ^ 
is trivially satisfied. If all the v in (fTUI) are given by then the previous proof applies. 
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4 The edge localizing potentials 



In this section we define the edge locahzing potential and show that its sum over the 
triangles is non-negative. The edge localizing potential is only nonzero for relatively flat 
triangles. Such triangles can have negative excess, but in the Delaunay triangulation they 
will be adjacent to a triangle with positive excess. The edge localizing potential exploits 
this fact. 

We will refer to the point equidistant to the three vertices of a triangle as the "cen- 
ter" of the triangle. Note that it need not lie inside the triangle. The edge potential 
e{xQ,Xi,X2,ri,r2) will depend on the signed distance from the center of the triangle to 
the edge opposite vertex 0. We let d{xo,Xi,X2) denote this signed distance. We define it 
to be positive if the center lies on the same side of the edge as vertex 0, and negative if 
they lie on opposite sides. We give the formula for this signed distance in appendix 1X1 
We note that the signed distance for an edge is negative when the angle opposite the edge 
is obtuse. 

Throughout this section we will consider two triangles that share an edge. The length 
of this common edge is Xq, and ri, r2 are the radii of the two discs centered at its endpoints. 
Xi,X2 are the lengths of the other two edges in one triangle and the lengths of 

the other two edges in the other triangle. 

Now consider (i(xo, Xi, X2) and (i(xo, x'^, Xg), the signed distances from the centers of 
the two triangles to their common edge. We claim that if the two triangles come from a 
Delaunay decomposition, then 



This inequality is obviously not true for arbitrary triangles that share an edge. It says 
that if the signed distance from the center of a triangle to an edge is negative, then the 
signed distance from the center of the other triangle that shares this edge to the edge 
must be positive and greater in absolute value. Given this inequality we can take any 
function f{d,xo,ri,r2) which is an increasing and odd function of d, and let 



Then §^ implies (fTT]). 

To prove ((221) ^^e figure |21 We have drawn the two triangles so their common edge 
is vertical, and we have shown the bisector of this common edge with a dashed line. The 
two centers are each equidistant from the endpoints of this common edge. So both centers 
lie on the dashed line. Inequahty ((221) is equivalent to C being to the right of C (or equal 
to C). By definition the three vertices of a triangle are equidistant to the center of the 
triangle. The Delaunay decomposition has the property that no other vertex of a triangle 
is closer to the center than these three vertices. In particular, C is closer to P than to P', 
and C is closer to P' than to P. It follows that C is to the right of C (or equal to it). 



d{xo, Xi,X2) + d{xo, x[, X2) > 



(32) 



e(xo,Xi,X2,ri,r2) = /(rf(xo, Xi, X2), Xq, n, rs) 



(33) 
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Figure 2: C and C are the centers of the left and right triangles respectively. In this 
example the signed distances from C and C to the common edge are negative and positive, 
respectively. 



A fairly simple function f{d,xo,ri,r2) will suffice for our purposes. Recall that each 
of ri and r2 is either r or 1. We define 



As always, f{d, xq, 1, r) = f{d, xq, r, 1). Finally, we let f{d, xq, 1, 1) = 0. 

5 Local proof the localizing potential works 

We now turn to the proof of ((Tj). We will use the computer to prove this inequality for 
most triangles. However, for triangles that appear in the densest packing shown in figure 
HJ equality holds in ((7j). Thus for triangles close to those that appear in this densest 
packing inequality ((Tj) will be close to an equality. So we can only use the computer to 
prove this inequality for triangles which do not lie too close to a triangle in the densest 
packing. In this section we will prove for the triangles which are close to a densest 
packing triangle. 




(34) 



Next we define 




(35) 
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There are four triangles that appear in this densest configuration. In all of these 
triangles the three discs are tangent to one another. As we will see this implies that the 
excess E{T) has a local minimum at each of these triangles. The localizing potential V{T) 
also has a local minimum at these triangles. We must show that when we perturb one of 
these triangles the increase in E{T) is greater than the increase in V{T). 

Throughout this section we will obtain bounds on quantities by bounding their partial 
derivatives with respect to Xo,Xi and X2. If f{xo,Xi,X2) is a function of the three edge 
lengths, and we have 

df 

Ci < -Q^{xo,xi,X2) < di (36) 

throughout some neighborhood of (xq, Xi, X2), then 

2 2 

CiAxi < f{xo, xi, X2) - f{xo, xi, X2) < ^ diAxi (37) 

i=0 1=0 

where Axi = Xi — Xi. 

We first consider triangles that have discs of radius r at all three vertices and which 
satisfy 

2r<Xi<2r + e, i = 0,1,2 (38) 

where e > will be determined later. Since the three discs all have radius r, D is 
independent of the Xi. Letting 

AE = E{xo, Xi,X2, r, r, r) — E{2r, 2r, 2r, r, r, r), (39) 

then 

2 2 
AE = A{SA-D)>SY, «rrr(e) Ax, = (5a„.,(e) ^ Axi (40) 

i=0 1=0 

where 

Q-rrrle) = mm 7- — (41) 

OXi 

with the min over the above set of triangles. Note that this min is independent of i. 
Inequality (ffTj) in the appendix gives a lower bound 

We must compare the above with the increase in the localizing potential V{T). Let 
A(f)i denote the change in the angles (pi corresponding to changing each xj by Axj. Then 
the change in V{T) is bounded by 

2 2 

AV{T) < m ^ \A(j)i\ < m ^ K„{e)Axj (42) 

i=0 j=0 

where 

^rrr(e) = max ^ 1—^1 (43) 

i=0 ^ 
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and the max is over the set of triangles. &rrr(^) is independent of i, and we denote it by 
brrri^)- This proves that the excess is greater than the locahzing potential for all triangles 
satisfying (j38|) provided m < rrirrri^) where 

mrrr{e) = 6^^ (44) 

For triangles that have discs of radius 1 at all three vertices the estimates are very 
similar. We consider the set of triangles which satisfy 

2<Xi<2 + e, i = 0,1,2 (45) 

The constraint we obtain on m is now m < miii(e) where 

Oiii(e) 

aiii(e) and &iii(e) are defined by analogy to arrr(e) and 6rrr(e)- 

Next we consider triangles with one disc of radius 1 and two of radius r. We take the 
large disc to be at vertex and consider the set of triangles given by 

2r<Xo<2r + e, 1 + r < Xj < 1 + r + e, z = 0,1,2 (47) 

In this case D is not independent of the Xj. Using the constraint 0o + 0i + 02 = ^r, we 
can write D as 



1 1 — 

D = -i(i)o + r^(pi + r^h) = + ^—^0 (4J 



Thus 



2 ^ „2 



= A(M -D)>J2 5<.(e) Ax, - ^^A0o (49) 

1=0 



where a\^^ is the min of 4^ over the set of triangles. 



dxi 

Let 



Cirr(e) = mm — 



di„(e) = max — — 

OXi 



(50) 

where the min and max are over the triangles satisfying (jTfjl . Note that there are no 
absolute values on the partial derivatives in the above. The calculations in the appendix 
show that on the set of triangles is positive for i = and negative for i = 1,2. We 
now have 



ae>y: 



i=0 



1 — r 



2 
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Ax, (51) 



We bound the increase in the locahzing potential V{T) as before. The constraints we 
obtain on m are m < m\^j.{e) with 



m 



Irr \ 



(52) 



Finally we consider triangles with one disc of radius r and two of radius 1. We take 
the small disc to be at vertex and consider the set of triangles given by 



2<xo<2 + e, l+r<Xi<l + r + e, 



0,1,2 



We now have 



D 



1 



(r^0o + 01 + 02) 



IT 



r^-l 



Thus 

2 -I 2 ^ 

i=0 i=0 

and so we obtain the constraints m < m*^]^(e) with 



6a. 



rll 



e + 



-Til 



(53) 
(54) 

(55) 



Kii(6) + i^4n(e) 



rll^ 



(56) 



^rll^ 



Using the results of appendix |Bl we can evaluate these bounds on m. We will take 
e = 0.001. For comparison we also show the values for e = 0. We find 

rrirrriO) = 0.135463, m,.„(0.001) = 0.134576 
mni(O) = 0.455814, miii(O.OOl) = 0.454183 
m?„(0) = 0.232960, m?^^(0.001) = 0.231100 
mi^^(O) = 0.179205, m*i^^(0.001) = 0.178067 
m°ii(0) = 0.264015, m°ii(0.001) = 0.262815 
m;ii(0) = 0.308628, m^^i (0.001) = 0.306788 

(57) 

where z is 1 or 2. Recall that we take m = 0.12, so these bounds are all easily met. In 
fact, with m = 0.12 this local proof works for e as large as 0.018. 



6 Global proof the localizing potential works 

We now prove inequality (|7j) for the rest of the triangles, i.e., all the triangles that are 
not close to those that appear in the densest configuration in figure [T] We use interval 
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arithmetic for this proof. This is necessary for two reasons. First, it is needed (as in 
all computer-assisted proofs involving real numbers) to deal with the fact that computer 
calculations involving real numbers are not exact. Second, even if the computer could 
perform exact calculations we could not prove ((Zj) by checking one triangle at a time since 
there are infinitely many triangles. Instead we must work with sets of triangles. Given 
intervals for the three edge lengths Xq, Xi, x^ of the triangle, we consider the set of triangles 
whose edge lengths belong to the respective intervals. It is straightforward to use interval 
arithmetic to then compute intervals for quantities such as the angles 0O) 0i) 02- Interval 
arithmetic is designed so that the interval computed for 0j means that for a triangle whose 
edge lengths belong to their respective intervals, the value of 0i must belong to its interval. 
We note that the code for computing intervals for quantitites such as 0, is the same as 
the code for computing these quantities when we have a single real value for each edge 
length. We need only change the data type of the variables involved from "double" to 
"interval" and define versions of the basic arithmetic operations (addition, multiplication, 
inverse cosine, etc.) for interval variables. The localizing potential involves real numbers 
that must be represented by intervals, and so it too is handled using interval arithmetic. 

We now fix choices of tq, ri, r2 and ask what triangles must be considered in the proof 
of (|7j). Since discs cannot overlap, we have the following lower bounds on the length of 
the edges of the triangle. 

Xq > ri + r2 
^2 > ^0 + 

(58) 

Recall that the center of a triangle is the point equidistant to the three vertices, and 
the circumradius is the radius of the circle centered at this point which contains all three 
vertices. In the Delaunay decomposition no disc has a center closer to the triangle's center 
than the discs at its vertices. So if the circumradius is greater than 1 + r, then we can 
add a disc of radius r to the packing by putting its center at the center of the triangle. 
If we assume that the packing is saturated, then this implies that every triangle has a 
circumradius of at most 1 + r. The circumradius is greater than half the length of any 
side of the triangle. So we have 

a;i<2(l + r), z = 0,l,2 (59) 

At each stage of the computer proof we have a list of parallelepipeds Jq x Ji x I2. Such 
a parallelepiped represents the set of triangles with Xi E li. At each stage of the computer 
proof the union of the parallelepipeds in the list is the set of triangles for which we must 
still prove Initially, there is one parallelpiped in the list with 

/o = h + r2,2(l + r)], Ji = [ro + r2,2(l + r)], = [r^ + r,,2{l + r)] (60) 
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At each stage we take a parallelepiped from the list and split it into two parallelepipeds. 
We compute the interval for E{T) — F{T) for each of the two. If it contains no negative 
values for a particular parallelepiped, then we know that ((Zj) is true for all triangles in 
this parallelepiped, and we can discard it. If the interval contains some negative values, 
we must add this parallelpiped to the list. Note that if the parallelepiped is contained 
in one of the parallelepipeds (gS)), (g?!), or then we know E{T) - F{T) > 

on this parallelepiped, and so it can be discarded. In appendix ^ we show that if the 
packing is saturated then in the Delaunay triangulation the area of every triangle is at 
least 2r^/(l + r). So for a parallelepiped, if the interval we compute for the area of the 
triangles is entirely greater than this lower bound, then we can discard the parallelepiped. 
The process ends when the list of parallelepipeds for which (jTj) has not been proved is 
empty. 

We compute our interval bounds in the simplest, crudest way. For example, we com- 
pute the interval bound on the angle 0j by simply using interval arithmetic in the equation 
for 0j in terms of the edge lengths, eq. (jHT|) . One could compute a better interval bound 
by explicitly maximizing and minimizing 0j over the parallelepiped by considering the 
derivatives of (pi with respect to the edge lengths Xj. Because of our rather crude bounds, 
the parallelepipeds must be split a large number of times. The computer considers a total 
of about 26 million parallelepipeds. This takes about 30 minutes on a laptop. 

A Formulae 

The vertices are labelled 0,1,2. The length of the edge opposite vertex i is xi. The radius 
of the disc at vertex i is r^. Each can only be r or 1. The angle at vertex i is 0j. It is 
given by 

cos(0i) = [x] + xl- x\)l {2xjXk) (61) 
where {i,j, k} = {0, 1, 2}. The area of the triangle is 

Recall that by the center of a triangle we mean the center of the circle which contains 
the three vertices of the triangle. We will derive a formula for the signed distance from 
the center to each edge and a formula for the circumradius. We take the vertices of the 
triangle to be 

Po = (0,0), Pi = (xo,0) P2 = {a,b) (63) 

where + 6^ = X2- Let P be the center of the triangle. The center must lie on the line 
X = Xo/2, so P = {xo/2,do), where do is the signed distance from the center to the edge 
opposite vertex 0. 

The circumradius of the triangle, R, satisfies 

i?= IIP -Poll = IIP- P-^ll (64) 
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This implies 

i?2 = ^ ^2 _ _ 2,^/2)2 + (6 - doY (65) 
Using + b'^ = x^, this gives 

axo + 2bdo = x\ (66) 



Using 



I 

a = X2 cos((/)i) = a;2 ° ^ ^ (67) 

2x0X2 



and 

\hxo = A (68) 

we find 

do = ^0(^1+ ^2 -a:g) ^gg^ 

8A 

Using ()65p. the above formula for (ig yields a formula for the circumradius. 

R=^J^ (70) 
For a saturated packing we can assume 

R<l + r (71) 

Otherwise we could add a disc of radius r with center at the center of the triangle. Using 
the above formula for R this implies a lower bound on the area A. 

A > (72) 
-4(l + r) ^ ' 

Using the trivial bound Xj > 2r, this implies 

A>-— (73) 
1 + r 

B Calculations for triangles with excess near zero 

We consider triangles with xq G [x, x + e], Xi, X2 G [y, y + e]. 
We have 

-LD-/i. — j^X qX ~t~ ZjX qX<2 ~r ^X ^X <2 Xq X-|^ X2 i J 

Given an i = 0, 1 or 2, we let j and k be the other two integers in {0, 1, 2}. Then 

dA 

8A— = x,{x^^+xl-x^) (75) 
17 



We assume that e is small enough that the set of triangles we are considering only contains 
acute triangles. So x| + — xf > 0. Then ^ > for each i, and so A is increasing in 
each Xi. So A is bounded above by its value for the triangle with sides of length xq = x + e, 
xi = X2 = y + e. This area is 

An^a. = \{x + e)V% + e)2-(x + e)2 (76) 



Thus we obtain the lower bounds 



and for i = 1,2 



From ()6ip we have 



OA ^ x{2y'-{x + ey) ^^^^ 

dXQ 8Amax 



dA y(x^ + y^-(y + eY) 



(sin(0o))"^ xj + xl-xl ^g^^ 



dxi 2x1X2 

= 2xlx, ^''^ 

Note that ^ is always positive and the other two partials are negative for acute triangles. 
So sin(0o) is bounded below by its value at xq = x, xi = X2 = y + e. Let 5*0 be the value of 
sin((/)o) at this point. And sin(0o) is bounded above by its value at xq = x + e, xi = X2 = y. 
Let 5*0 be the value of sin((/)o) at this point. Then we have 

X ^d^^x + e 



(y + eyS'o - dxo " y'S, 







for 2 = 1,2. 
Similarly 



{y + e)^ + (x + - ^ ^ ^ y^ + x^ - {y + 

2y^So - dx, - 2{y + eyS', ^ ^ 



1 fsin(0i))"' ^ (84) 



dXi XqX2 
■*! / ■ / A Xq -\- X^ — X2 /or-\ 

- (sm(0i)) (85) 



dxo 2xqX2 



- (sm(0i)) (86) 



dx2 2x1x0 
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For these derivatives we only need bounds on their absolute values. 



-I < (87) 



dxi xySi 
d(j), ^ {x + ef + {y + ef-y^ 

dct^, 2{y + ef - x' 

' ' - 2y'^xSi ^ ' 

with equal to the value of sin(0i) when xi = y, xq = x + e and X2 = y + e. 

C Equation for r 

In this appendix we show that the radius r for the compact packing is a root of an eighth 
degree polynomial. The angles satisfy 

2a- (3 = 120 (90) 

So 

cos(q; - 13/2) = cos(60) = 1/2 (91) 

Hence 

cos(q;) cos(/?/2) + sin(Q;) sin(/?/2) = 1/2 (92) 

This imphes 

sin2(Q;) siT?{l3/2) = (1/2 - cos(q;) cos{l3/2)f (93) 

which is equivalent to 

1 - cos^ {af - cos^{(5/2) = J - cos(q;) cos(/?/2) (94) 

We have 

, , 2(l + r)2-22 r2 + 2r-l 

= 2(1 + = ITTTF 

and 

sin(/?/2) = ^ (96) 

cos(/3/2) = ^^^^ (97) 
Substituting these into the above equation, some algebra gives 

- 8r^ - 44r^ - 232r^ - 482r^ - 24r^ + 388r^ - 120r + 9 = (98) 
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By considering the packing itself, the conditions that various pairs of discs are tangent 
can be used to show that r also satisfies 

(7 + AVsy + (20 + uVsy + (6 + AVsy + (-20 - 4V3)r + 3 = (99) 

To compute the density of the packing, consider the unit cell shown in figure H The 
quadrilateral has sides of length 

/ = 1 + r + VSr + Vl + 2r (100) 

and angles of tt/2 and 2tc/3. So its area is -\/3/^/2. It covers a total of 3 large discs and 3 
small discs. So the density is 

6 = ^^^(^ + ^') ^ 0.911627478 (101) 
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